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Fano resonanceslT], the signature of multi-path inter- 
ference between a continuum and discrete resonance, have 
been exploited in opticsiS & for a variety of appUcations 
including biosensingiT], optical switching|5] and wave- 
length conversion|6]. In cavity-QED|7], photons stored in 
a cavity leak into a continuum of radiation modes along 
with the atomic spontaneous emission. Measurements of 
this radiation can, in principle, yield information regard- 
ing the strength of atom-cavity interactionslisi [§]. Utiliz- 
ing nano-assembly techniques to integrate the constituent 
components of a soUd-state cavity-QED system, here we 
realize a platform for studying interference phenomena 
of an emitter coupled to a microcavity and its radiation 
mode environment. A quantum model of the system is 
presented, from which the coherent coupling rate between 
cavity and emitter is estimated. It is envisioned that this 
nano-assembly approach may be appUed to the integra- 
tion of more complex cavity-QED geometries, where the 
optical coupling and entanglement of multiple quantum 
emitters can be reaUzed. 

The cavity QED system studied in this Letter (Fig. [T]i 
consists of a diamond nanocrystal containing optically ac- 
tive nitrogen-vacancy (NV) centers that is coupled to a di- 
electric microdisk cavity. Optical access to the nanocrystal- 
microdisk system is provided by both an optical fiber taper 
waveguidel 10| and a high numerical aperture (NA=0.55) ob- 
jective lens. In bulk diamond, NV centers exhibit atom-like 
optical properties, and are promising candidates for appli- 
cations in quantum information processing, includingsingle 
photon generation! 1 1]> coherent population trapping|l2lj and 
optical readout and manipulation of single nuclear spinsllO, 
[l4l [Tsll . NV centers hosted in diamond nanocrystalsflg] 
are amenable to nanomanipulation[17J and integration with 
nanophotonic structures lll8[|l9ll . The fiber taper waveguide is 
a versatile tool in this regard; in addition to serving its usual 
role as a probe for the optical fields of the microcavity, in this 
work it also functions as a collection optic and a positioning 
tool for diamond nanocrystals. This provides several bene- 
fits, notably pre-selection of nanocrystals with optimal spec- 
tral properties from a nominally lower quality ensemble 1 2^, 
and the controlled positioning of any number of nanocrystals 



on photonic structures. 

Light radiated from a coupled emitter (hereafter referred to 
as a "dipole") and cavity system is usually described by two, 
distinct, source channels. The cavity radiation channel con- 
sists of the localized quasi-mode of the cavity and the radia- 
tion modes which it leaks into, whereas the dipole radiation 
channel consists of all other radiation modes directly coupled 
to the dipole emitter. In the system studied here, the radiation 
mode emission from the cavity mode and the dipole overlap, 
as illustrated for a generic Fabry-Perot cavity in Fig. [He), re- 
sulting in interference between the radiated field from each 
channel. A quantitative understanding of the interference ef- 
fects can be obtained through a quantum mechanical model of 
the dipole-cavity system. For simplicity we consider a single 
mode of the microdisk which leaks into radiation modes with 
energy decay rate 2k. The NV center optical transition is mod- 
eled by a single dipole transition, with excited state energy 
decay rate due to spontaneous emission and pure dephasing 
rate Jp. The coherent-coupling rate between the dipole transi- 
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FIG. 1: (a) Schematic of the nanocrystal fiber-to-microcavity posi- 
tioning technique, (b) SEM image of a diamond nanocrystal posi- 
tioned on the edge of a Si02 microdisk. Illustration of the interfer- 
ence between direct dipole and indirect dipole-cavity radiation, (c) 
for a generic Fabry-Perot single-sided cavity and (d) for a microdisk 
cavity with nanocrystal-scattering radiation loss. 
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tion and an excitation of the microcavity mode is g. The elec- 
tric field radiated into the far-field by the dipole-cavity system 
can then be written as. 



E(r,f) =e,/(r)e-"P''Vy;a_(f)+e,(r)e-*V2Kfl(f)+h.c., 

(1) 

where a is the microcavity field operator and d_ is the polar- 
ization operator of the dipole transition. In general, the field 
profiles ed{r) and ec(r), of radiation from the dipole and the 
microcavity, respectively, need not be orthogonal. The phases 
(^d,c associated with the field operators a and d- depend on the 
system and measurement geometry. In the generic example 
drawn in Fig. [Uc), the relative phase between the direct and 
indirect emission is determined by the additional path length 
followed by the cavity emission, and by the phase imparted 
by the cavity output coupling. For the nanocrystal-microdisk 
system (Fig.[TJd)), the path length difference can be zero. 

The time evolution of the operators d_ and a can be cal- 
culated from the system density matrix equation of motion, 
which depends on the dipole-cavity Hamiltonian and on Lind- 
blad operators representing the cavity and dipole decoherence 
processes[9] (see Appendix). In the "room temperature" limit, 
Yp ^ K,Ys,g, Aco, and the detected optical spectra into a given 
collection optic is: 
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where Aco — (Oc — (£) is the cavity-dipole emission detun- 
ing, £d,c — (eo(r)|ef/x (r))r describes the overlap between the 
mode e„ (r) of the collection optic and that of the dipole (cav- 
ity) radiation. Note that the relative amplitude of the cavity 
emission scales with — Ig^/Kjs, the bad-cavity Purcell fac- 
tor Also note that in the room temperature limit, the effect of 
phonon-assisted emission on the dipole-cavity dynamics can 
be included as a dominant contribution to jp. 

The microdisk cavities studied in this work are 20/jm in 
diameter, 400 nm thick, and are formed by thermal oxi- 
dization of pre-patterned Si microdisks ll2lll . The resulting 
Si02 microdisks (refractive index «si09 ~ 1 -45) support high- 
Q whispering-gallery-modes (WGMs) across the visible and 
near-infrared spectrum, with measured values as high as Q = 
8 X 10^ in the 830-855 nm wavelength band. Of particular in- 
terest are the transverse-magnetic-like (TM-like) disk modes, 
which are polarized primarily normal to the planar microdisk 
resulting in significant electric field intensity at and above the 
disk surface (Fig. |4tb,c)). For a wavelength of 637 nm, near 
the emission wavelength of the negatively charged NV transi- 
tion (NV^), the effective mode volume of the lowest radial- 
order ip = 1) TM-like WGM (TMp=i) is calculated to be 
Veff = 82(A-/nsiOT)^- The maximum ratio between the elec- 
tric field energy density at the disk surface to that at the point 
of peak electric field energy density (lying within the disk) is 
l/r| = 0.23. For the NV^ transition of a nanocrystal placed 
at this optimal surface location, and with dipole orientation 
aligned normal to the disk surface, this translates into a co- 
herent coupling rate between the TMp=i mode and the NV^ 



FIG. 2: Scanning confocal microscope (SCM) images (details in 
[l2]) of a mesa before (a) and after (b) a nanocrystal has been picked 
up with a fiber taper. Optical images of (c) fiber taper and attached 
nanocrystal aligned with microdisk edge, and (d) after nanocrystal 
placement. Nanocrystal imaging is aided by a white light source 
coupled into the fiber taper and the microdisk, (e) SCM image of 
microdisk after nanocrystal placement. 



dipole of go/^n « 0.64 GHz. This estimate is based upon 
a total excited state spontaneous emission rate of 'Y\\/2k Ki 12 
MHz ll22ll . The situation, however, is complicated for the NV^ 
transition due to electron-phonon interactions which result in 
signigicant coupling to several phonon sidebands. Of partic- 
ular interest is coupling to the NV^ zero phonon line (ZPL), 
which for the 3-5% fraction of spontaneous emission that is 
emitted into the ZPL i23ll24ll . yields a reduced coherent cou- 
pHng rate to the TMp=i of gzpi/2Jt « 0.13 GHz. 

The nanopositioning scheme used to assemble individual 
diamond nanocrystals on the surface of the microdisks is il- 
lustrated in Fig. [n^)- A nanoscale taper formed from single 
mode silica fiber (Nufern 630-HP) is used both as a collection 
optic and as a means to pick-and-place diamond nanocrys- 
tals. The diamond nanocrystals (Diamond Innovation-NAT, 
200nm median diameter) are initially sparsely dispersed on 
a clean silicon sample patterned with elevated mesa struc- 
tures. The diamond coated silicon sample, the microdisk sam- 
ple, and the optical fiber taper are mounted in an enclosed 
box with a dry nitrogen purge. After identifying a nanocrys- 
tal of interesting !2ff], the fiber taper is contacted on top 
of a nanocrystal using high resolution (50 nm) stages. The 
taper is then raised vertically away from the mesa surface. 
Due to surface interactions, which are not yet fully under- 
stood, the nanocrystal attaches to the fiber taper and is lifted 
from the silicon mesa surface with high yield (see Figure|2fa- 
b)). Transferring the nanocrystal to a microcavity follows the 
"pick-up" process in reverse (Fig. |2c-d)). In order to de- 
tach the nanocrystal from the fiber taper it was found nec- 
essary to slightly move the microdisk (1 — 10 /jm in-plane 
steps) when in contact with the nanocrystal so as to rub it off 
of the taper. Images of the microdisk after placement of the 
nanocrystal (« 200 nm diameter) are shown in Figs. [Hb) and 
Fig. |2je). Using this technique nanocrystals can typically be 
placed within a taper diameter (< 500 nm) of the disk edge. In 
this case, the nanocrystal has been placed approximately 150 
nm from the disk edge. To facilitate widefield optical imaging 
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FIG. 3: (a) Emission from a diamond nanocrystal attached to a fiber 
taper, collected through the high-NA lens in the far-field and through 
the fiber taper in the near-field. The fiber taper data was scaled by a 
factor of 1.6 to take measured fiber insertion loss into account. Peaks 
near X = 690 nm are due to fluorescence in the Ge doped fiber core, 
(b) Measured emission when the taper interacts with a microdisk, for 
varying nanocrystal position relative to the microdisk and spectrom- 
eter, as indicated by the illustrations. 



of the nanocrystal during positioning, a relatively large parti- 
cle was chosen in this instance; smaller nanocrystals can also 
be manipulated and positioned appropriately with the aid of 
confocal microscopy to identify the nanocrystal position. 

At each step in this pick-and-place process we use the fiber 
taper, in combination with widefield imaging optics, to study 
the nanocrystal emission. Figure [3|a) shows the photolumi- 
nescence (PL) spectra after the first step in which the dia- 
mond nanocrystal is attached to the fiber taper. Emission is 
collected both from a high-NA objective lens at normal inci- 
dence to the sample as well as through the fiber taper. Op- 
tical excitation is performed using a frequency-doubled YAG 
laser (X w 532 nm). The pump laser beam is sent through 
the collection-lens and focused down to a 1.5 /jm spot, with 
a resulting pump intensity of ~ 10 kW/cm^. Emission from 
the NV^ ZPL is visible at A, « 637nm, superimposed on the 
large phonon sideband characteristic of diamond NV centers. 
The measured efficiency of the fiber taper collection relative 
to the lens collection is ^ 40%. Theoretically, using smaller 
diameter fiber tapers 300nm) and nanocrystals 50nm), 
it should be possible to reach an absolute taper collection effi- 
ciency of > 10% 125,26]. 

Figure[3jb) shows measured PL spectra during the next step 
in the process in which the nanocrystal is brought into the 
near-field of the microdisk using the fiber taper In this step 
the fiber taper waveguide is aligned and evanescently coupled 



to the microdisk with the nanocrystal positioned (i) far from 
the microdisk, on the nearside of the spectrometer input, (ii) 
nearly touching the microdisk in the taper-microdisk coupling 
region, and (iii) far from the microdisk, on the far side of the 
spectrometer input. In (i) the microdisk has no significant ef- 
fect on the measured NV emission spectrum as only light di- 
rectly emitted into the fiber reaches the spectrometer In (iii) 
the microdisk behaves as a simple drop filter on the emission 
radiated into the fiber, resulting in Lorentzian dips at each of 
the cavity resonances. The spectrum in (ii) is more complex, 
with various Fano lineshapes appearing in the spectrum due to 
interference between the cavity and taper spontaneous emis- 
sion channels (more about this below). 

In the final step of the nanocrystal-microdisk assembly pro- 
cess, the nanocrystal is transferred onto the surface of the mi- 
crodisk. The strong interaction between the nanocrystal and 
microdisk WGMs is evidenced by the effect of the nanocrys- 
tal on the WGM spectral properties. Prior to placement of 
the nanocrystal, the microdisk in Fig. [Ttb) supports a pair 
of degenerate TEp=i traveling wave WGM resonances with 
e = 3.4 X 10^ in the 852 nm wavelength band (blue curves 
of Fig. Ua)). After placement of the nanocrystal near the 
disk edge, the WGM resonance splits into an asymmetric 
pair of standing wave resonances (red curves of Fig. Ufa)) 
formed from backscattering by the nanocrystal [27, 28]. The 
standing wave cavity modes spatially lock to the position of 
the sub-wavelength nanocrystal, with the anti-node (node) 
of the lower (higher) frequency resonance aligned with the 
nanocrystal. As a result, the g-factor of the long wavelength 
resonance is degraded to Q — 1.7 x 10^ due to scattering by 
the nanocrystal, while the Q-f^ctor of the shorter wavelength 
resonance remains relatively close to its unperturbed value. 
The measured mode-splitting and scattering loss induced by 
the nanocrystal is in good correspondence with a perturba- 
tive analysis (see Appendix), depending upon similar cavity 
field properties as the dipole-cavity coherent coupling rate g 

iHH]. 

Optical excitation of the assembled nanocrystal- 
microcavity system shows the striking influence of the 
microdisk cavity, and the resulting interference between 
the different channels for nanocrystal emission. PL spectra 
collected in the far-field using the high-NA lens and the fiber 
taper are shown in Fig. |5] The fiber taper is evanescently 
coupled to the microdisk on the edge diametrically opposite 
from the nanocrystal, and only collects light emitted into the 
cavity modes of the microdisk. In both collection geometries 
the envelope of the collected spectrum follows the broad 
NV^ emission characteristic of Fig. [3jb); we estimate that 
the nanocrystal under study here contains less than 5 NV 
centers. The fiber taper PL spectrum consists primarily 
of regularly spaced Lorentzian peaks corresponding to the 
microdisk WGM resonances. The lens-collected PL, a 
high-resolution spectrum of which is shown in Fig. |3b), is 
instead punctuated by Fano-like resonances superimposed 
upon the broad background spectrum of the NV^ transition. 

For a nanocrystal positioned on the microdisk top sur- 
face, 150 nm from the disk edge, FEM simulations show 
that in the 600 nm wavelength band the higher radial order 
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FIG. 4: (a) Microdisk mode lineshapes measured by monitoring tiie 
fiber taper transmission spectrum with the taper positioned in the mi- 
crodisk near-field, before (blue) and after (red) nanocrystal place- 
ment. A 850nm band tunable diode laser source (New Focus Veloc- 
ity) was used to measure the taper transmission; a 630nm tunable 
source was not readily available at the time of the measurements. 
The dashed lines are fits | 21] to the resonant modes. Cross-sectional 
view of the 630 nm wavelength band (b) TMp=i and (c) TMp=3 
whispering-gallery mode profiles near the disk edge as calculated 
using the finite-element-method (FEM). The white outline indicates 
the periphery of the microdisk, with the shape of the disk profile esti- 
mated from SEM images. Only the dominant, vertical component of 
the electric field is plotted for clarity. The position of the nanocrystal, 
as placed on the disk, is indicated by the red diamond. 
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FIG. 5: (a) PL from the diamond nanocrystal placed on the mi- 
crodisk in Fig.[nb), collected at room temperature through the near- 
field fiber taper (green) and far-field lens (blue), (b) High resolution 
(6A. = 20 pm) PL spectrum of the lens-collected emission around 
X = 680 nm. The red curve is a fit using 5iens from the text, general- 
ized to include multiple decoupled cavity modes. 



(p = 3,4) WGMs of the microdisk are most strongly cou- 
pled to the nanocrystal (Fig. |4tb,c)). A perturbative anal- 
ysis of the scattering loss induced by the nanocrystal (see 
Appendix) also indicates that the Q values of the p — 1-4 
radial order modes, save the TMp=4, should be limited by 
nanocrystal scattering. From careful comparison of the mea- 
sured emission spectra with the calculated radiation-limited 
and nanocrystal-scattering-limited Q values, the various fam- 
ilies of cavity modes in the PL from the nanocrystal can be 
identified, as indicated in Fig. |5] Although visible in the 
high-resolution lens-collected PL spectrum (Fig. |3b)), the 
TEp=3 4 modes are not faithfully reproduced due to the lim- 
ited resolution of the spectrometer. Note that doublet splitting 
is not expected for the these modes, owing to their lower-Q 
and large mode volume (see Appendix.) Absent from both 
the lens and taper collected emission spectra are the highest 
Q, p — 1,2 radial-order modes, a result of the low spectrome- 
ter resolution which washes out narrow spectral features, and 
the tighter confinement of these modes inside the microdisk 
which weakens the coupling to the fiber taper waveguide. 

To understand the measured PL spectra from the two collec- 
tion geometries further, we now turn to the model described by 
eq. (|2]). For collection through the fiber taper, we observe no 
interference due to negligible overlap between the fiber taper 
mode and the radiation modes of the nanocrystal (ej/ec ^ 0). 
The resulting emission spectrum is then dominated by the sec- 
ond term in eq. (|2| resulting in Lorentzian peaks centered at 
each of the microdisk cavity resonance wavelengths as seen 
in the measured spectrum of Fig.|5] For the far-field lens col- 
lection, e^/ w £( when the cavity radiation is dominated by 
scattering from the nanocrystal, as is estimated for all but the 
TMp=4 modes. For a sub-wavelength nanocrystal in which the 



origin of the dipole emission and the scattered cavity emission 
nearly coincide, one can show that (|)^ — = 7t/2 (see Ap- 
pendix), resulting in a lens-collected power spectral density, 



5lens(tt') °= — 




(3) 



Using the above relation for ^lens. a fit to the Fano reso- 
nances of the high resolution lens-collected spectrum is per- 
formed (Fig. |5lb)). For the TMp=3 modes, which are both 
spectrometer-resolved and scattering-limited, the fit yields an 
estimate of k/2k = 15 GHz and Fo 0.20. In the case of the 
spectrometer-resolved TMp=4 modes we find k/2k = 73 GHz 
and Ff, w 0.020; however, the reduced kFo product for this 
mode family is consistent with a radiation-limited Q-factor 
for which < £j. 

Most directly, F„ represents a ratio of the local density of 
states of the high-Q cavity modes relative to that of the ra- 
diation modes, as seen by an NV center embedded in the 
nanocrystal. Fitting the data gives the ratio between and 
js, from which one can obtain g by substituting the appro- 
priate Yj. Of particular interest is the discrete ZPL transition 
of the NV^ state whose line width can approach its radiation- 
lifetime-limited value of 12 MHz at cryogenic temperatures 
m. Substituting (y,. = Yzpi)/2:i: ~ 0.5 MHz, k/2k = 15 GHz 
and Fa — 0.2 for the TMp=3 resonances yields gzpi/27t w 28 
MHz. This value is roughly half that estimated from the sim- 
ulated TMp=3 mode overlap with the nanocrystal (68 MHz), 
and within a factor of 5 of the estimated maximum value 
for this microdisk (130 MHz). Factors contributing to the 
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smaller observed ^^pi include imperfect mode-matching of 
the direct dipole emission and the nanocrystal-scattered cavity 
field, local field effects stemming from polarization of the di- 
amond host, and the NV position(s) within the relatively large 
nanocrystal, for which the presence of multiple NV centers 
with random dipole orientations is likely to result in a lower 
average g^^i. 

These measurements provide an initial demonstration of de- 
terministic placement of single diamond nanocrystals on opti- 
cal microcavities. Using this technique in future experiments 
to selectively manipulate smaller nanocrystals containing sin- 
gle, narrow linewidth NV centers, it should be possible to 
fabricate a coupled NV-cavity system with Purcell enhanced 
emission exceeding unity. Such a system would be an impor- 
tant step towards the efficient optical readout and control of 
the electronic and nuclear states of an NV center, and the re- 
alization of NV center based devices for quantum information 
processing f?^, such as a quantum repeater More impor- 
tantly, this method may be utilized to integrate more complex 
diamond nanocrystal cavity QED systems, forming the basis 
of a chip-based quantum network. 



Methods 



Veff^ 



/«^(r)|E,(r)pjV 
max[«2(r)|E,(r)|2] 



(Al) 



where n(r) is the refractive index profile of the cavity and 
Et(r) is the electric field of the cavity mode. For the 
whispering-gallery-modes of the microdisk, the quoted effec- 
tive mode volumes, unless otherwise stated, are for standing- 
wave resonances, which are a factor of two smaller than their 
traveling-wave counterparts. We use a factor r|(ro) to relate 
the peak of the electric field energy density in the cavity to 
the electric field energy density at a particular position To 
(^{^o) > !)■ The per photon electric field amplitude at po- 
sition To then relates to the effective mode volume as. 



photon (I'o ) 



2e„n2(r„)r|(r„)yeff' 



(A2) 



from which the coherent coupling rate can be estimated as 
g = I • Ec photon (fo) I /fi- The magnitude of the dipole moment 
of the electronic transition, can be determined from the 
spontaneous emission rate in bulk material, 



Microdisk fabrication The microdisk optical cavities are Si02 disks, 20 /jm 
in diameter, approximately 400nm thick near the disk edge, and supported 
by a central Si02 pillar ~ 5 fjm in diameter. They were fabricated by ther- 
mally oxidizing Si microdisks in an oxygen purged furnace at a temperature 
of 1050°C. The template Si microdisks were fabricated from a siHcon-on- 
insulator (SOI) wafer with a 2 /jm buried-oxide-layer and a p-doped 217 nm 
Si device layer of resistivity 14-20 Clem. The Si microdisk processing was as 
follows. Electron beam Hthography was used to define the microdisk pattern 
in a polymer electron beam resist (ZEP 520A). To improve circularity and 
remove surface roughness, the patterned resist was reflowed at 160°C [2^. 
A SF6/C4F8 inductively coupled reactive ion etch transfers the resist pattern 
into the top Si layer An HE wet etch partially removes the underlying Si02 
layer, resulting in Si microdisks supported by Si02 posts. 
Photoluminescence measurements Diamond nanocrystals containing NVs 
are excited optically using a frequency-doubled YAG laser (X « 532 nm). A 
50X (NA = 0.55) ultralong working distance microscope objective lens fo- 
cuses the excitation beam to a 1.5 fjm spot. Typical excitation power is 250 
/.iW. A dichroic mirror on the backside of the objective separates NV photo- 
luminescence from excitation power Photoluminescence is passed through 
a long wavelength pass filter (cut-off ~ 540 nm), and then directed with a 
flipper mirror to either a spectrometer (resolution ~ 20pm), or imaging optics 
and a CCD camera. 
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APPENDIX A: EFFECTIVE MODE VOLUME AND 
COHERENT COUPLING RATE 

In this work we use the commonly defined effective mode 
volume of a microcavity mode in terms of its peak electric 
field energy density; 



2 37t2feoC^Y| 



«ncCO-' 



(A3) 



APPENDIX B: MODE-COUPLING FROM A 
SUB-WAVELENGTH SCATTERER 

The modal coupling between clockwise (cw) and coun- 
terclockwise (ccw) traveling wave modes in a whispering- 
gallery-mode microcavity has been observed experimentally 
and explaine d by many other authors, including those of Refs. 
ll2ll I27I [30I I31I [32I1 . Here, we present a simple analysis of 
this coupling. Maxwell's wave equation for the vector electric 
field in a microdisk structure is 



V^E^ 



("0 



0, 



(Bl) 



where /jq is the permeability of free space, e** is the dielectric 
function for the ideal (perfectly circular) microdisk and 5e is 
the dielectric perturbation that is the source of mode coupling 
between the cw and ccw modes. Assuming a harmonic time 
dependence, the modes of the ideal microdisk structure are 
E^(r,f ) = E^(r)exp(/C0yf), and are solutions of eq. dSTT ) with 
5e = 0. Solutions to eq. ( IBlb with 5e 7^ (i.e., modes of the 
perturbed structure) are written as a sum of the unperturbed 
mode basis 



Eir,t)=e-"'»''£ajm%r). 



(B2) 



Plugging this equation into eq. ( IBll i. and utilizing mode or- 
thogonality, we arrive at the coupled mode equations 
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coo/5£(E»(r))*.E»(rMr 
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(B3) 
(B4) 



where Aco^t = coa — coo- In deriving eq. ( IB3I ) we have assumed 
that the mode ampUtudes change slowly relative to coo, and 
that jcOy — cool <C COo- We have also ignored small "self term" 
corrections to the eigenfrequencies. Reference |21] presents 
a functional form for p in situations involving small surface 
roughness perturbation. Under weak scattering conditions an 
assumption is made that only each pair of localized, degener- 
ate cw and ccw WGMs of azimuthal mode number ±m (with 
dominant polarization (TE or TM) and radial mode number 
p) are coupled by the disk perturbation 5e. The coupled mode 
equations for these traveling wave modes then read as 



dt 

dciccw 

dt 



-iA(£)a,„{t) + i\^\e'^accw{t) 
- iAmoccw it) + i\^\e''^acw{t), 



(B5) 



with P = I P I e'^ given by the integral in cylindrical coordinates 



. _ coo JiJe+'^"'^8e{<^,p,z)d<^) (E°,(p,z))'p^p^z 
P 2 /eO|EO„,(r)Pdr ' ^ 

where we have used the fact that E°^,(r) = E°„,(p,z)e+""'^ and 
Etcvi,(r) = (E[lj,,(r))*. For a sub-wavelength nanocrystal scat- 
terer, the dielectric perturbation in eq. (IB6I 1 can be approxi- 
mated as 



5£ = 8o(«nc-l)8'''(r-rnc)Kc 



150 nm from the disk edge, the ratio of the mode-splitting to 
the center frequency of the TEp=i WGMs is measured to be 
8.2 X 10^^. This corresponds well with the theoretical value 
predicted by eq. ( IB8I 1 of 2.2 x 10^^ for a simulated traveling- 
wave effective mode volume Vm-, eff = 86(A,„/«si02)^ ™d r] = 
0.024 evaluated at the center of the nanocrystal (100 nm above 
the disk surface). 



APPENDIX C: SURFACE-SCATTERING FROM A 
SUB-WAVELENGTH SCATTERER 

Optical losses in whispering-gallery-mode resonators are 
often limited by refractive index perturbations, 5e, on the cav- 
ity surface. These index perturbations are sourced approxi- 
mately by the unperturbed field solutions to create polariza- 
tion currents. 



/coSeE. 



(CI) 



In analogy with microwave electronics, the polarization 
currents drive new electromagnetic fields which radiate into 
space. Optical losses due to the perturbations can be calcu- 



lated from the far field solutions setup by J 
vector potential sourced by J is given bvll33 



The far field 



A™.(r) = ^ 



,-ikQr 



j 3{r')e-'''°''' dr' , (C2) 



(B7) 



where is the wave vector in the surrounding air, and we 
have made the simplifying assumption that the microcavity 
does not significantly modify the dipole radiation from its free 
space radiation pattern. 

For a point-like nanocrystal scatterer with perturbation 
given by eq. ( IB7I ). and for source field corresponding to a 
standing-wave WGM of the microdisk, the far field Poynting 
vector is given by 



where Eq is the free space permittivity, Wnc is the refractive 
index of the diamond nanocrystal 2.4) and Vnc is the phys- 
ical volume of the nanocrystal. The mode-splitting in angular 
frequency is 2|P|, and is proportional to the center frequency 
COq. The normalized mode-splitting, in terms of the traveling- 
wave mode effective volume (note this is twice as large as the 
standing-wave mode volume), can then be written as 



1 _2|p| _ (nge-l)Kc 
gp COo ri(rnc)Vr,y,eff ' 



(B8) 



where r|(rnc) is the correction factor taking into account the 
reduced electric field strength at the position of the nanocrys- 
tal. 

We have performed measurements of the scattering effects 
on the TEp=i WGMs in the 850 nm wavelength band of the 10 
/jm radius, Si02 microdisks studied in this work. For the 200 
nm diameter nanocrystal placed on the microdisk top surface. 



Sw = r- — \rxKrad\ 



(C3) 



)V„2eo|E,(rnc)|' |f xe(rn 



0)1^ 



327t2 r2 

where E^ (rnc) is the cavity mode electric field and e(rnc) is the 
unit vector representing the polarization of the cavity mode 
electric field at the position of the nanocrystal scatterer The 
total (cycle-averaged) power radiated, Pmci, far from the disk 
can be found by integrating the outward intensity over a large 
sphere. 



Prad 



(S-f) r^d^ 
&kl (n 



2 -l)V2eo|E,(r„ 



327t2 



(C4) 

■J |f xe(rnc))pfl'n. 
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TABLE I; Calculated mode parameters in the X ~ 600 nm band. TABLE IL Calculated mode parameters in the A. ~ 850 nm band. 



Mode 


Gmd 






rinc 




m 


"eff 






(^o/"SiO, 


)^ 










TEp=i 


10'^ 


69 


0.057 


0.013 


3.6 X lO'* 


125 


1.27 


TM^=l 


10" 


82 


0.061 


0.021 


2.7 X 10"* 


122 


1.24 


TEp=2 


lOlO 


97 


0.11 


0.026 


2.5 X 10"* 


121 


1.23 




108 


106 


0.21 


0.069 


1.0 X 10"* 


118 


1.20 


TE^=3 


1.7 X 10^ 


103 


0.12 


0.028 


2.5 X 10"* 


115 


1.17 


TM^=3 


8.6 X 10^ 


106 


0.24 


0.079 


9.0 X 10^ 


113 


1.15 


TEp=4 


1.7 X 10^ 


109 


0.11 


0.026 


2.9 X 10"* 


111 


1.13 


TMp=4 


1.7 X 10"* 


111 


0.23 


0.077 


1.0 X 10"* 


109 


1.11 



Mode 


Qi-dd 




^.s' ^nc Qss 


"eff "2 












TEp=i 


4.5 X 10** 


43 


0.073 0.024 1.4 X 10^ 


1.27 93 


TMp=i 


2.2 X lO*" 


51 


0.086 0.040 8.4 x 10^ 


1.22 87 



APPENDIX E: NANOCRYSTAL-CAVITY SPECTRUM 
CALCULATION 

The measured spectrum for the nanocrystal-microdisk sys- 
tem is given by 



For the sphere centered about the nanocrystal scatterer, 
/ |f X e{rnc))f' dQ. = 87t/3. The quality factor of a cavity is 
given by g = OiUc/Prad, where Uc ^ \J 8''(r) lEc pt/r is the 
cycle-averaged stored energy in the cavity. The cavity energy 
can also be related to the effective mode volume through eq. 
(lAlb . Combining all of these relations, we can rewrite eq. 
(IC4b as a scattering quahty factor 



e.v = 



3X^11 (rnc)yeff 
47t2(«2^-l)2y2 



(C5) 



Estimates of the scattering quality factor for the various 
WGMs studied in this work are presented in Tables HI and HH 
where the nanocrystal has a measured diameter of 200 nm and 
the correction factor r|(rne) is evaluated at the center of the 
nanocrystal (approximately 100 nm above the disk surface). 



APPENDIX D: FINITE-ELEMENT-METHOD 
SIMULATIONS OF SIO2 MICRODISK MODES 

In Table |T] and we present the results of finite-element- 
method (FEM) simulations of the 10 /jm radius, Si02 mi- 
crodisks used in this study at resonant wavelengths in the 650 
nm and 850 nm wavelength bands, respectively. The effec- 
tive index of each cavity mode is calculated from the approx- 
imate relation, «eff ~ mXo/lnRo, where m is the azimuthal 
mode number of the WGM and Rg = 10 /jm is the physical 
radius of the microdisk. The correction factors, r], and r|„t, 
correspond to the electric field energy density evaluated at the 
radial position of the nanocrystal 150 nm from the disk 
edge) and vertical position at the surface of the microdisk and 
at the center of the nanocrystal (100 nm above the disk sur- 
face), respectively. The surface-scattering-limited Q-factor is 
estimated from eq. (IC5l l using r[„c. 



5lens(C0) =Re 



dt I dt e' 




,ico(f-r') 



E:(r,f)-E„(r,f') 



(El) 

where E„(r,f) is the component of the radiated nanocrystal- 
microdisk field, E(r,f), collected by the microscope objec- 
tive or fiber taper collection optic. E(r,f) is given as a func- 
tion of the system variables a{t) and d_(r) by eq. (1) in the 
manuscript text. To predict the spectrum of a spontaneous 
emission event, we assume that at f = 0, the dipole is pre- 
pared in its excited state and the microdisk contains no pho- 
tons, and follow the techniques given in Refs. |l9l|34|] to cal- 
culate 5iens(K>) from (a(f)) and (a_(f)), as a function of the 
coupling and loss parameters of the system, K, yj , y^, and de- 
tuning Aco between the dipole and microdisk resonances. An 
outline of this method is given below. 

The system Hamiltonian, which describes the dipole-cavity 
interaction in absence of decoherence, is given by 



(E2) 



where coj and co^ are the dipole transition and cavity mode 
frequencies, d_ = and a are the dipole and cavity excita- 
tion annihilation operators, respectively, and g is the coherent 
dipole-cavity coupling rate defined in [A] Coupling between 
the system variables and the measured radiation field is taken 
into account using the quantum master equation for the system 
density matrix p: 



U. + Lp)p 



where i-a.cp are Lindblad operators given by: 



- (2a_pd+ - 6+d-P - p6+d_ ) 



L^p — K(2flpfl^ — fl'ap — pfl^a) 
Y 



(E3) 

(E4) 
(E5) 
(E6) 



with = d+d_ — d_d+. L4 (Lc) describes decay of the 
dipole (cavity) excitation into radiation modes, at energy de- 
cay rate y( (2k). Lp represents decoherence of the dipole ex- 
citation, in the form of pure non-radiative dephasing at rate 
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Equations of motion for the classical expectation values of 
d and d_ can found from the quantum master equation, using 

the relation 



dt 

d . 



-(/co,. + k) (fl)+^(d_) 



(E7) 
(E8) 



Note that in absence of an external driving field, when a dipole 
is initially prepared in its excited state in a cavity with zero 
photons, (fl d;) = — (fl) for all time f. These equations of mo- 
tion can be used calculate the two-time correlation functions 
required to determine the optical spectrum defined by eq. (lEll ). 
This is accomplished using the quantum regression theorem, 
which, from eqs. ( IE7b and (IE8I ). allows us to write, for t' > f , 



— {aHt)a{t'))^-{mc + K){a\t)a{t')) 



(E9) 



dt 



y(fl+(f)a-(f')) = + f+yp)(fl'(f)d_(f')) (ElO) 



-8 {a 



\t)a{t')) 



dt 



-{d+{t)a{t')) = -{i(ac + K){^+{tW)) + g{d+{t)&^{t')) 

(Ell) 



dt 



-(d+(f)d-(f')) = -(/co^ + |+yp)(6+(f)a_(f')) (E12) 

-g{&+{t)a{t')). 



Writing these equations, as well as analogous equations of 
motion for the special case f = f', in the frequency domain, 
a closed set of algebraic equations for the various contribu- 
tions to 5iens (h)) is obtained. Note that we have not assumed 
that the spontaneous emission process is stationary. In the 
room-temperature limit, Yp >• [ | CO — cOf/.f | , ^, K, Ys] , we find 



(fl-(f)d_(0> 

{a'{t)a{t')) 



= 



1 



Q^((0)=F[(d+(f)d_(0)] = 
Q,(co)=F[(a+(f)«(0)] = 



1 



1 



YpYv + 



(E13) 
(E14) 
(E15) 
(E16) 



where F[g(r,r')] = £ dt £ dt' e'"^^'-''^ g{t ,t') and Aco = co- 
cOc, and we set co^/ = cOf. From eqs. (|E13l - |ET6] l, eq. (IE II ). and 
eq. (1) in the text, we arrive at the expression for ^(co) given 
by eq. (2) in the text. 



APPENDIX F: DIPOLE-CAVITY RADIATION PHASE LAG 

Here we derive the phase difference between the radiation 
directly emitted from a dipole source (e.g., a diamond NV 
center) and the scattered radiation from a cavity resonance 
driven by the same dipole source. In particular we consider the 
situation explored in the main text in which the cavity scatter- 
ing site is subwavelength and coincident with the dipole site 
(e.g., a diamond nanocrystal). We will show that the phase dif- 
ference, 8(|) = (|)^/ — is simply due to the 7t/2 phase shift be- 
tween driving source and local cavity field when driven on res- 
onance. We begin with Maxwell's equations for a microcavity 
defined by £t(r), a dielectric nanocrystal perturbation Ae„(r) 
which couples the microcavity mode to radiation modes, and 
a dipole source term Js(r,f): 
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VxVxE(r,f) + ^e,(r)E(r,f): 



^Ae„(r)E(r,f)-|;J,(r,f) 
(Fl) 



We have assumed that the microcavity is non magnetic, and 
that the vacuum dielectric and magnetic permittivities are 
unity. 

To derive temporal equations of motion, we expand E in 
terms of the cavity and radiation modes. 



E(r , f ) = c(f )e, (r)e-'"'' + ^ O (f )e, (r)e- 

j 



(F2) 



where et.y(r) are solutions to Maxwell's equations in absence 
of the nanocrystal and external source. 



V X V X ecj(r) - co^ 8c(r)ecj(r) = 0, 



(F3) 



e;(r) 


= 0, 


(F4) 


e,(r) 


= hi, 


(F5) 


e.;(r) 


= 1. 


(F6) 



and (£)j,c are real. It follows that ec,/(r) are orthogonal, and we 
choose to normalize the fields as follows: 



('1e.(r)|;) 



Note that in this mode expansion, et(r) represents a truly 
bound mode of the cavity: in absence of the nanocrystal or 
other perturbation, excitations of this mode have an infinite 
lifetime. 

Substituting eq. (IF2b into eq. (IFll l. and applying eqs. ( IF3I 
- IF6b . for a harmonic driving term Ji(r,f) = Jo(r)e^'"'', the 
following coupled mode equations can be obtained: 



W r. <-J J\ J r. 



dt 



dt 
where. 



o(0 = 



2co,' ^ 



-1 ft),-o)iU 



^cj^K%^{c\Aen\j), 



(F7) 
, (F8) 



(F9) 



9 



(c,;'|Jo) 



(FIO) 



In deriving eqs. (IF7b -( IF8l l. we have ignored "self-terms" 
which result in corrections to the eigenfrequencies, and have 
assumed that the mode amplitudes vary slowly compared to 
the optical frequencies. We have also assumed that the field 
ec(r) is predominantly confined within the cavity, and that 
fields (r) are predominantly excluded from the cavity. Upon 
integrating eq. ( IF8I 1 directly, and substituting the result into eq. 
(IF7b . we arrive at a differential equation for c(f). Assuming 
that the radiation modes form a continuum with slowly vary- 
ing density of states, p„, and coupling coefficients, = Kcj, 
and that coupling from the radiation modes into the cavity 
is negligible, a simplified equation of motion for c(f) is ob- 
tained: 



^C{t) = (-K + /(C0,-(0,))c(f) + ^i,, 



(Fll) 



where c{t) = e'("'""')'c(f ), and 2k = 7ip„co^ | p/2 is the total 
energy decay rate of the cavity mode into the radiation mode 
continuum, due to scattering from the nanocrystal. Averaging 
over long times, only radiation into the mode whose frequency 
is resonant with the source frequency will be appreciable; the 
equation of motion for the mode amplitude, = rj=s, is 



dt 



co,.K*. 



2 4 i{(£)c - COj 



(F12) 



where fj{t) = r,(f)e''"' and we have used the steady 
state solution for c{t) from eq. (IFl II) . 



The Fano nature of radiation from the nanocrystal-cavity 
system can seen in the righthand side eq. (IF121 i: the first term 
represents emission from the source "directly" into the radia- 
tion modes, while the second term represents "indirect" emis- 
sion from the source into the radiation modes via the cavity 
mode. Note that eq. (IF12b has the same form as eq. (2) in 
the text, the primary difference lying in the definitions of the 
various coupling coefficients. Here, the relative phase of the 
coupling coefficients is well defined by Maxwell's equation, 
and can easily be determined in the case of a subwavelength 
scatter such as a nanocrystal. 

Assuming that Ja(r) = j„8'^'(r — r„) and Ae„(r) = 
5e„5'^^ (r — r^), the coupling coefficients are given by, 

K, = e:(r„) •e,(r„)5e„ - |e:.(r„) •e,(r<,)5e„| , 

(F13) 

s,. = -e:(r„) • j„ = |e*(r„) • j„| , (F14) 

= -e:(r„) • j„ = |e:(r,) • j„| , (F15) 

and eq. (IF12b . with the relative phases included explicitly, is 



\Sc\ i(OcK-\ 
4 /((Oc-co.s) 



(F16) 



The above equation indicates that for the system of interest, 
when the cavity and source are on resonance, "direct" and 
"indirect" contributions to the radiation field are k/2 out of 
phase. 
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